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ABSTRACT 

We present a study of the evolution of the fraction of radio-loud active galactic nuclei (AGN) 
as a function of their host stellar mass. We make use of two samples of radio galaxies; one in 
the local universe, 0.01 < z ^ 0.3, using a combined SDSS-NVSS sample and one at higher 
redshifts, 0.5 < z ^ 2, constructed from the VLA-COSMOS_DEEP Radio Survey at 1.4 GHz 
and a K^-selected catalogue of the COSMOS/UltraVISTA field. We observe an increase of 
more than an order of magnitude in the fraction of lower mass galaxies (M* < Mq) 

which host Radio-Loud AGN with radio powers ^ W Hz“^at z ^ 1 — 2 while 

the radio-loud fraction for higher mass galaxies (M* > 10^^ Mq) remains the same. We 
argue that this increase is driven largely by the increase in cold or radiative mode accretion 
with increasing cold gas supply at earlier epochs. The increasing population of low mass 
Radio-Loud AGN can thus explain the upturn in the Radio Luminosity Eunction (RLE) at 
high redshift which is important for understanding the impact of AGN feedback in galaxy 
evolution. 

Key words: galaxies; active - radio continuum; galaxies - galaxies; evolution - galaxies; jets 
- galaxies; luminosity function, mass function - accretion, accretion discs 


1 INTRODUCTION 


During recent years it has become increasingly apparent that the 
radio jets of radio-loud Active Galactic Nuclei (RL AGN or radio- 
AGN) play a crucial role in the process of galaxy formation and 


evolution via ‘AGN feedback’ (see e.g. |Best et al. 

o 

o 

To 

o 

o 

-o 

[Bower et al.[[2006| [Croton et al. [[20061 [Fabian et al. 

[2006k This 


feedback can occur because the enormous energy output of the 
AGN can be injected into the surrounding medium, possibly also 
the fuel source of the AGN, via ionizing radiation and/or relativis¬ 
tic jets, thereby providing enough energy to affect star formation in 
the host galaxy. Radio galaxies have been shown to comprise two 
different populations: high and low excitation ( [Best et ar] | 2005[ 
|Tasse et al.|2008[ |Hickox et al.|2009| l, each of which may have a 
separate and different effect of feedback, the exact nature and evo¬ 
lution of which is still debated. The first population of RL AGN 
is associated with the classic optical ‘quasars’. These sources ra¬ 
diate across the electromagentic spectrum and are consistent with 
the unified models of quasars where emission is obscured at some 
wavelengths when the source is viewed edge-on (e.g. |Barthel|1989| 
|Antonucci|1993[[Urry & Padovani|1995^ . In this ‘high-excitation’, 
‘cold mode’ or ‘radiative mode’, accretion is postulated to occur 
via an accretion disc in a radiatively efficient manner (e.g. |Shakura| 
|& Sunyaev|1973^ . These high excitation radio galaxies (HERGs) 

are typically hosted by lower mass, bluer galaxies in less dense en¬ 
vironments (e.g.|Tasse et al.|2008[|Janssen et al.|2012^. The second 
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mode of radio activity was first noted by their lack of strong high- 
excitation narrow-line optical excitation expected from the ‘quasar’ 


mode ([Hine & Longair|1979[[Laing et al.|1994[ [Jackson & 

Rawl-[ 

[ings[ 1997k Moreover they show no evidence of mid-infrared emis- 

sion from dusty tori IWhysong & Antonucci|2004||0gle et al. 

20061 

and no evidence of accretion-related X-ray emission ([Hare 

castle[ 


|et al.|2006{ [Evans et al.|2006) . [Hardcastle et alT]j2007| first sug- 


gested that this mode, known as the ‘low-excitation’, ‘radio mode’, 
‘hot mode’ or ‘jet mode’ occurs when hot gas is accreted directly 
onto the supermassive black hole in a radiatively inefficient manner 
via advection dominated accretion flows (ADAEs, e.g.|Narayan &| 
[YT||1995k|Best et ^j2005t showed that these low excitation ra¬ 
dio galaxies (LERGs) are hosted by fundamentally different galax¬ 
ies: higher mass, redder and occurring in more dense environments. 
This mode in particular provides a direct feedback connection be¬ 
tween the AGN and its hot gas fuel supply in the manner of work 
done by the expanding radio lobes on the hot intra-cluster gas. Eor 
a more detailed review of the HERG versus LERG dichotomy see 
[Heckman & BestH2014[l and references therein. 


In order to understand the relative significance of the differ¬ 
ent types of radio-AGN feedback we need to understand the cos¬ 
mic evolution of radio sources in detail. It has been known for 
several decades that the comoving number density of powerful ra¬ 
dio sources is two to three orders of magnitude greater at a red- 


shift of two to three compared to the local Universe (e.g. Schmidt 


1968 

Sandage[[1972| Osmer[[1982| Peacock 

1985 

[Schmidt et al. 

1988 

Dunlop & Peacock|1990| Rigby et al. 

2011 

. Similarly, the 


space density of optically selected quasars (QSOs) peaks at red- 
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shift 2 < 2 < 3 (e.g. |Boyl e et al.|198 8{ |Hewett et al.|19^|Warr^ 
|et al.|[T9M| >. It is also well known that, within the local universe 
(z < 0.3), the fraction of galaxies which host a radio source, i.e. 
the radio-loud fraction, is a very s teep function of h ost galaxy stel¬ 
lar mass (/radio-loud « 


Best et al. 


20051, 

. i i - J 


increasing to 

> 30 per cent at stellar masses above 5 x ICP'Mq for radio lu¬ 
minosities > 10^^ WHz“^. This pervasiveness of radio-loudness 
among high mass galaxies in the local Universe, combined with the 
dramatic increase in density of radio-loud sources at earlier times 
suggests that there must be an increase in the prevalence of radio 
activity among galaxies of lower mass. To test this, we investigate 
the fraction of radio-loud sources out to redshift ~ 2 as a function 
of their host stellar mass. 

This paper is structured as follows. In Section we describe 
the consh'uction of the Radio-Loud AGN samples and the cata¬ 
logues from which they are selected. In Section the luminosity 
functions for these samples are constructed and binned by host 
galaxy stellar mass. The radio-loud fraction as a function of host 
stellar mass is also determined. We discuss the results and their 
implications in Section]^ and conclude in SectionThroughout 
this paper, we use the following cosmological parameters: Hq = 
70kms“^ Mpc“^, = 0.3 and Da = 0.7. The spectral index, 
a, is defined as Sv oc and unless otherwise specified, we adopt 
a default value of 0.8. 


2 RADIO-LOUD AGN SAMPLES 

To investigate the evolution of the Radio Luminosity Function 
(RLF) for different stellar-masses over a similar luminosity range 
several samples with good ancillary and derived data are needed, 
to provide both the required statistics at low redshift and the sen¬ 
sitivity at high redshift. In this work we combine one already ex¬ 
isting matched radio-optical dataset in the local universe using the 
SDSS-NVSS sample, described in Section|2T] and one at redshifts 
0.5 < z < 2, which we have constructed using the VLA-COSMOS 
DEEP Radio Survey at 1.4 GHz and a Ks-selected catalogue of the 
COSMOSAJltraVISTA field, which is described in Section [T2l 


2.1 Local SDSS Sample 

Eor the local radio source sample we use the catalogue compiled by 
|Best & Heckman|j2012|, which was constructed by cross-matching 
optical galaxies from the seventh data release (DR7; |Abazajia n| 
|et al.|200^ of the Sloan Digital Sky Survey (SDSS) spectroscopic 
sample with radio sources in the NRAO Very Large Array (VLA) 
Sky Survey (NVSS; [Condon et al.|1998^ and the Eaint Images of 
the Radio Sky at Twenty centimetres (FIRST; [Becker et al.|1995^ . 
The parent optical sample consists of all galaxies in the value- 
added spectroscopic catalogues (VASC) created by the Max Plank 
Institute for Astrophysics and Johns Hopkins University (MPA- 
JHU) group (see [Brinchmann et al.||20()4[ available at http: // 
WWW .mpa-garching. mpg. de/SDSS/I. The cross-matching 
was done for all radio sources with flux densities > 5 mjy, which 
corresponds to radio luminosities of P-^ 4 GHz ~ WHz“^at 
redshift 2 = 0.1 and P-^ 4 GHz ~ 10^"^ WHz“^at redshift 2 = 0.3. 
The combined radio-optical area covered is 2.17 str ( [Best & Heck- 
maiJ20]^^ Of the 927 522 galaxies in the VASC, [Best & Heck- 
man| j201^ selected a magnitude-limited sample of 18 286 ra¬ 
dio sources, which they showed to be 95 per cent complete and 
99 per cent reliable ([Best et al.|2005^. The sample was restricted to 


the ‘main galaxy sample ’ jStrauss et al.|2002) , comprising galax¬ 
ies within the magnitude range 14.5 ^ r < 17.7mag and the 
redshift range 0.01 < 2 ^ 0.3. This local radio-optical sample 
consists of 9168 radio sources. We note that, being based on the 
SDSS main galaxy sample, this local matched radio-optical sam¬ 
ple excludes both radio-loud quasars and broad-line radio galaxies. 
However, [Best et aL[j2014) show that this is only problematic at 
radio powers above P-^ 4 GHz ~ WHz“^. Since our LFs do 
not probe those high powers we make no correction for this bias. 
Moreover, we know that the dominant population of radio sources 
in this sample are not quasars. 

Properties of the host galaxies are taken from the VASCs 
which, for each source, provide several basic measured parame¬ 
ters from the imaging data such as magnitudes, colours and sizes 
( [York et al.|[200Q| ), as well as derived properties including, most 
importantly for this work, the stellar mass M* ( Kauffmann et al.j 
|2003| l. For their matched radio sample, [Best & HeckmanH2012 1 
also separated the sources into star-forming galaxies and RL AGN 
(7302 sources), which are further sub-divided into high-excitation 
(HERG) and low-excitation (LERG) sources, based on their optical 
photometric and spectroscopic parameters. 

The left panels in Eig.[T]show the radio power and host stel¬ 
lar mass as a function of redshift for the radio sources in the 
SDSS-NVSS sample. All the radio sources are plotted in black, and 
sources in the restricted sub-sample within the redshift bin which 
defines our local sample are plotted in yellow. 


2.2 Distant VLA-COSMOS Sample 

The VLA-COSMOS Large Project covered the 2 deg^ of the COS¬ 
MOS Eield at 1.4 GHz with observations by the VLA in the A con¬ 
figuration. This survey, extensively described in jSchinnerer et aL[ 
( |2004[[2007| >, provides continuum radio observations with a reso¬ 
lution (half-power beam width) of 1.5 x 1.4 arcsec and a mean Icr 
sensitivity of about 10.5 /rjy in the innermost 1 deg^ region and of 
about 15 /iJy in the outer parts. Prom the VLA-COSMOS catalogue 
we selected sources above 50/rJy, which corresponds to sources 
> 5(7 over 50 per cent of the survey area. 

The optical data for this sample comes from the Ka-selected 
catalogue of the COSMOS/UltraVISTA field ( jMuzzin et al.|2013| l 
which contains PSF-matched photometry in 30 photometric bands 
covering the wavelength range 0.15pm —>■ 24pm. The entire re¬ 
gion overlaps with the radio survey so the combined area is that of 
the COSMOS/Ultra VISTA data, 1.62 deg^. Pollowing the recom¬ 
mended criteria of |Muzzin et al.| ( [2013t , we selected a ‘clean’ sam¬ 
ple from the Ks catalogue of sources with flags: ‘star’ 7 ^ 1, ‘K_flag’ 
^ 4, ‘contamination’ ^ 1, and ‘nan_contam’^ 3 (these flags relate 
to stars and saturated sources and the quality of the photometry for 
nearby sources). Pinally we selected sources brighter than the 90 
per cent completeness limit of Ks ^ 23.4 mag. The redshifts pro¬ 
vided are determined from SED fitting to the broadband photom¬ 
etry using the EAZY code ( jBrammer et al.|2008^ and for source 
2 < 1.5 they quote an rms error of {Sz/{1 + 2 )) = 0.013 and 
a catastrophic outlier fraction of 1.56 per cent. While the quoted 
outlier fraction is very low, we do expect that many of the outliers 
will be quasars (i.e. HERGs). It is well known that photomePic 
redshifts determined for quasars are less reliable than those ob¬ 
tained for galaxies (e.g. RichaMs^talj^OH |Babbedge et al.|2004t 


[Mobasher et al.|20^|Polletta et al.|2007t. 

Within this sample we defined three ‘high-redshift’ sub¬ 
samples in the redshift ranges 0.5 < 2 ^ 1.0, 1.0 < 2 1.5 

and 1.5 < 2 ^ 2.0, chosen such that each contains ~ 300 — 700 
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Figure 1. (left) Full SDSS-NVSS Sample and the selected ‘local’ 0.01 < z 0.3 sample : radio power vs z and stellar mass, M», vs z. (right) Full VLA- 
COSMOS Sample and the three selected ‘high redshift’ 0.5 <2^1 and 1 < 2 : ^ 2 samples : radio power vs z and stellar mass, M*, vs z. The grey 
shaded regions show the regions in mass-space, redshift-space and radio power-space which demarcate the redshift samples, mass bins and power bins for the 
Luminosity Functions. The selected clean mass-complete radio samples are shown in yellow, orange and red symbols (these colours denote these redshift bins 
thought this paper). 


galaxies. The right panels of Fig.[^show the radio power and host 
stellar mass as a function of redshift for the radio sources in the 
VLA-COSMOS sample. Sources in the restricted ranges which de¬ 
fine our high redshift samples are plotted in orange and red, and all 
the remaining sources are shown in black. 

Note that at radio powers 4 GHz ~ 10^^ WHz“^, the sam¬ 
ples should consist almost entirely of RL AGN - we expect very 
little contamination from star-forming galaxies. We confirmed this 
using the spectroscopy-based AGN/SF separation in the local sam¬ 
ple ( [Best et al.|2005^ . Moreover, the st ar formation ra te for this ra¬ 
dio power is in excess of 25 Mq yr“^ 1 Condon 1992 1 so we expect 
AGN to continue to dominate even at higher redshifts. 


3 THE STELLAR-MASS DEPENDENT LUMINOSITY 
EUNCTION 

3.1 The Luminosity Eunction 

A standard technique for quantifying the rate of evolution of a pop¬ 
ulation of galaxies is to compare their luminosity functions (LFs) 
at two different epochs. In this section, we therefore determine the 
evolution of the Radio-Loud AGN population among host galaxies 
of different masses by comparing the luminosity functions of the 
SDSS-NVSS and VLA-COSMOS samples. Determining the radio 


LFs for the full samples first serves to confirm our sample selection 
and methods. 


These radio luminosity functions were calculated in the stan¬ 
dard way, using p = TjiljVi method 1 Schmidt| 1968 Condon 


|1989| >, where Vi is the volume in which a given source could be de 
tected. This volume is determined by both the minimum and max¬ 
imum distance at which a given source would be included in the 
sample as a result of the selection criteria: Vi = Umax — Utnin’ 
where Vfnax and are the volumes enclosed within the ob¬ 

served sky area out to the maximum and minimum distances re¬ 
spectively. The minimum accessible volume is a result of the opti¬ 
cal cut-off on the bright end (> 14.5 mag for SDSS-NVSS, while 
the VLA-COSMOS sample has no bright flux limit). The maxi¬ 
mum accessible volume is determined by the flux limits of both the 
optical « 17.77mag for SDSS-NVSS and < 23.4mag for VLA- 
COSMOS) and radio data (> 5 mjy for SDSS-NVSS and > 50 piy 
for the VLA-COSMOS sample). In practice, the maximum redshift 
is dominated by the radio flux limits. The normalization of the lu¬ 
minosity function requires knowledge of the precise intersection 
area of all input surveys. The sky area for the SDSS-NVSS sample 
is taken to be 2.17 sr jBest et aL]|2005| l. For the VLA-COSMOS 
sample, the sky area is taken to be 1.62 deg^ for all sources above 
100 ^Jy. For the faintest sources in the VLA-COSMOS sample, 
the area in which the sources could be detected is smaller due to 
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Table 1. Number of sources in each stellar mass-redshift bin. 


Mass bin 
log(M./M0) 

0.01 - 0.3 

Redshift bin 

0.5-1.0 1.0-1.5 

1.5 - 2.0 

10.00 - 10.75 

1021 

296 

231 

147 

10.75 - 11.00 

843 

153 

126 

90 

11.00 - 11.25 

1639 

127 

111 

74 

11.25 - 12.00 

5306 

51 

35 

26 

Total 

9006 

672 

518 

343 


the non-uniform rms noise level in the VLA-COSMOS mosaic. We 
therefore weight each source by the inverse of the area in which 
it can be detected, which also accounts for the varying detection 
area within a given luminosity bin. Uncertainties are calculated as 
the statistical Poissonian errors with a contribution of cosmic vari¬ 
ance where appropriate; in some luminosity/mass bins, these er¬ 
rors are so small that the error will be dominated by systematic 
effects. The area covered by VLA-COSMOS sample of 1.62 deg^ 
is small enough that the effects of cosmic variance are not negli¬ 
gible for very massive galaxies. The total area covered by COS- 
MOS/UltraVISTA is approximately 1.5 square degrees in one sin¬ 
gle field, making the effects of cosmic variance not negligible for 
very massive galaxies. The contribution of cosmic variance to the 
total error budget was estimated through the recipe of ( |Moster et al.| 
|2011| |. The average uncertainties due to this effect for the least mas¬ 
sive galaxies vary from 6 to 7 per cent and for the most massive 
galaxies vary from 10 to 12 per cent across the three highest red- 
shift bins. These values were added in quadrature to the Poissonian 
error of the LFs. 

We have compared our derived RLFs for radio-loud sources of 
all stellar masses with those in the literature. For the full samples, 
our local RLF agrees with that of |Best et al.| j2005^ and our RLF 
for the high redshift VLA-COSMOS sample agrees with that of 
ISmolcic et al.| ( [200^ . The agreement is not unexpected as we are 
using the same data, but it serves to validate out sample selections 
and method. 

We next split the samples into four bins of different host 
galaxy stellar masses: 10.0 < log(M*/M 0 ) ^ 10.75, 10.75 < 
log(M*/M 0 ) s; 11.0, 11.0 < log(M./M 0 ) ^ 11.25, and 
11.25 < \og{Mt/M q) sC 12.0. These bins are chosen to sam¬ 
ple the stellar masses well with similar numbers of sources in each 
bin (see Fig. [^. The RLFs derived in each of these stellar mass 
bins in the four redshift bins are plotted in Fig. Note that in 
the low redshift sample we exclude the points below P-^ ^ GHz ~ 
10^3 WHz“^for the highest stellar mass bin as these sources are 
only detectable out \o z < 0.1 and this sample is incomplete for 
high mass sources at these redshifts. Table [T] shows the number of 
sources in each stellar mass-redshift bin. From these LFs we see 
that in the local universe and for 4 GHz ~ WHz“^, the 
number density of the highest host stellar mass bin is the greatest, 
while in the higher redshift bins the number density within all the 
stellar mass bins is becoming increasingly similar. 

3.2 Space Density Evolution 

Now, to quantify the disproportional increase in radio activity 
among galaxies of different mass at higher redshifts, we investigate 
the relative comoving space density of radio-loud sources with re¬ 
spect to the local comoving space density. We do this by dividing 



log P [log W Hz 


Figure 2. The comoving space density of radio-loud sources. This is the 
RLF in redshift bins 0.01 < z ^ 0.3 (top panel), 0.5 < 2 ^ 1.0 (top- 
middle panel), 1.0 < 2 1.5 (bottom-middle panel) and 1.5 < 2 ^ 2.0 

(bottom panel) for the four host stellar mass bins (plotted in colour). 
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Table 2. Relative space density of radio-loud sources with respect to the lo¬ 
cal density as a function of redshift for all sources with radio powers greater 
than the cut-off Pj^ 4 GHz > lO^^WHz”^. 


Mass bin 
log(M./M 0 ) 

0.5 - 1.0 

Redshift bin 

1.0- 1.5 

1.5 - 2.0 

10.00 - 10.75 

8.9 ±3.1 

27.6 ± 6.9 

45 ± 11 

10.75 - 11.00 

9.5 ±2.2 

12.1 ± 2.1 

17.3 ± 2.6 

11.00 - 11.25 

5.7 ± 1.0 

6.45 ±0.86 

5.89 ±0.77 

11.25 - 12.00 

0.95 ±0.21 

0.68 ±0.14 



the stellar mass dependent RLFs in redshift bins 0.5 < 2 ^ 1.0, 
1.0 < 2 ^ 1.5 and 1.5 < 2 ^ 2.0 by the local RLF (0.01 < 2 ^ 
0.3). These relative RLFs are plotted in Fig. from which it is 
clear that there is a difference between the relative comoving space 
density of low mass galaxies hosting radio sources and that of high 
mass hosts. We derive the relative comoving space density of radio- 
loud sources with respect to the local comoving space density for 
all sources with radio powers greater than a cut-off luminosity of 
Pi 4 GP [2 ^ 10^^ WHz“^. This is plotted as a function of redshift 
in Fig.|^and the values listed in Table|^ 

At 2 ~ 1.5 — 2 the space density of the least massive galax¬ 
ies hosting Radio-Loud AGN above ^ GHz > IO^'^WHz-Ms 
45 ± 11 times greater than the local space density. Even the sources 
with stellar masses in the range 10.75 < log(M,/MQ) ^ 11.0 
are 17.3 ± 2.6 times more prevalent. We note that there is a slight 
decrement in the space density of the most massive galaxies hosting 
Radio-Loud AGN at both radio power cuts going out to 2 < 1.5. 
If we consider the slightly less powerful sources, P-^ ^ GHz ^ 
lO^S-S WHz“^, and only go out to redshift 2 < 1.5 (not plotted 
here), the effect is only really seen in the lowest mass bin where 
the increase in space density increases more rapidly with redshift: 
the space density of the least massive galaxies hosting Radio-Loud 
AGN is 40.8 ± 6.2 times greater than locally compared to the 
27.6 ± 6.9-fold increase for P^ 4 GHz ^ 10^"^ WHz~^sources in 
the 1.0 < 2 ^ 1.5 bin. 


4 THE RADIO-LOUD FRACTION 

Another way of looking at the increase in prevalence of AGN in 
lower mass hosts is to consider the fraction of sources which are 
radio-loud as a function of the host stellar mass in our four redshift 
bins. The mass-dependence of the radio-loud fraction can be an 
indicator of the accretion mode of the radio-AGN largely because 
of the different dependence of the fuelling source (hot vs. cold gas) 
on stellar mass jBest et al.|2006^ . The radio-loud fraction can be 
easily calculated by dividing the stellar mass function (SMF) for 
radio-loud sources, pnad (for sources above a given radio power 
limit), by the SMF for all galaxies, popt'- 

/radio-loud ~ PRadjPOpt- 


4.1 The Stellar Mass Function 


In order to calculate the radio-loud fraction we first derived the 
SMFs for all galaxies and for the radio source hosts by using the 
i/Vmax estimator 1 Schmidt|1968 as previously described for the 
LFs, which corrects for the fact that the samples are magnitude lim¬ 
ited. In order to construct the SMFs, the redshift-dependent limiting 
M* above which the magnitude-limited sample is complete needs 
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Figure 3. The relative comoving space density of radio-loud sources with 
respect to the local comoving space density. This is the RLF in redshift 
bins 0.5 < z 1.0 (top panel), 1.0 < 2 Si 1.5 (middle panel) and 
1.5 < 2 Si 2.0 (bottom panel) by the local, 0.01 < 2 Si 0.3, RLF for the 
four host stellar mass bins (plotted in colour). 


to be known. For the COSMOS/UltraVISTA sample we use the em¬ 
pirical 95 per cent completeness limit calculated by|Muzzin et al.| 
l |20T3t . 

The SMFs for radio-loud galaxies and all galaxies are plotted 

in Fig. 1^ for a radio-power cut-off of P 4 4 > 10 ^'*WHz 

We choose this limit because the highest redshift bin is only able to 
probe radio powers greater than this (see also Fig.[^. The SMFs we 
derive for our optical galaxy samples (for both the SDSS and COS¬ 
MOS samples) are consistent with those of |Muzzin et al.| j2013| l 
who use a more sophisticated maximum-likelihood analysis to de¬ 
rive the SMFs in several redshift bins. As expected ( [Best et aT] 
|2005[ >, pRad differs significantly from popt - the hosts of radio 
sources are biased towards more massive systems. Interestingly, 
while the comoving number density of all galaxies decreases with 
redshift, that of the radio source hosts with \og[Mt /Mq) < 11.0 
increases. This is consistent with the results of |Tasse et al.| j2008| l 
and shows that the radio-loud galaxy population evolves differently 
from the population of galaxies as a whole. 
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Figure 4. Relative space density of radio-loud sources with respect to the lo¬ 
cal density as a function of redshift for all sources with radio powers greater 
than the cut-off ^i.4GHz>10""WHz-i. 



logAf* [log Mg] 


Figure 5. SMFs for all galaxies, popt, (shaded lines) and radio-loud 
galaxies, pnad, (lines with points) for a radio-power cut-off ^ GHz ^ 
10^"^ WHz“^in the four redshift bins. 


4.2 The Radio-Loud Fraction 

Having computed the relevant SMFs we can calculate the radio- 
loud fraction as described above. Figure shows the radio-loud 
fraction for a radio-power cut-off of P-^ ^ GHz > lO^'^WHz"^ 
The radio-loud fraction clearly increases with redshift. Moreover, 
the slope of the mass dependence becomes shallower, showing 
that the fraction of lower mass galaxies hosting radio sources in¬ 
creases more with redshift than the fraction for high mass galax¬ 
ies. The slopes are shallower at higher redshifts, for ^ GHz ^ 
10^^ WHz“Agoing from /_r_l oc in the local sample, 

flattening to Jrl oc /rl oc and Jrl oc 

^i.oio.i jjj jjjg higher redshift bins. This is highlighted in the 
plot of the radio-loud fraction relative to the local redshift bin 
(fig- where the relative radio-loud fraction is up to two orders 
of magnitude greater at the low mass end. We note that the flatten¬ 
ing is quicker with redshift when the lower power P^ 4 GHz ^ 
10^®'® WHz“^sources are considered (in the first three redshift 



Figure 6. The fraction of galaxies hosting a radio source (radio-loud frac¬ 
tion) for a radio-power cut-off of ^ GHz > lO^^WHz ^ in the four 
redshift bins. The coloured lines show a linear fit over the stellar mass 
range 10 < log(M*/M 0 ) < 11.5. The slopes of these fits are 2.7 di 0.2, 
1.7di0.2, l.SdiO.l and 1.0 di 0.1 from the lowest to highest redshift bins. 



Figure 7. The fraction of galaxies hosting a radio source (radio-loud frac¬ 
tion) relative to the local fraction for a radio-power cut-off of ^ GHz ^ 
lO^'^ WHz“^in the three higher redshift bins. 


bins only, not plotted here): Jrl oc f^iL oc 

and/fli ocMi i±° i. 


5 INTERPRETATION 

In this section we aim to interpret and explain our results within 
the context of the HERG-LERG population dichotomy and their 
differential cosmic evolution and mass dependence. 

In the local Universe, the density of high mass Radio-Loud 
AGN is an order of magnitude higher than that of low mass Radio- 
Loud AGN at all radio powers (cf. top panel Fig. |^. In the more 
distant Universe, up to 2 < 2 we see a sharp increase in the num¬ 
ber density of Radio-Loud AGN hosted by lower mass galaxies, 
while the number density of high mass Radio-Loud AGN remains 
constant. 
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This large increase in the prevalence of radio activity among 
galaxies of lower mass at higher redshifts (cf. Fig. shows that 
it is the lower mass galaxies which are the cause of the upturn 
in the observed RTFs (e.g. |Dunlop & Peacock|1990[ |Rigby et al.| 
|2011| |. Moreover we suggest that this upturn is likely due to an in¬ 
creasing population of cold mode accretors at earlier epochs. From 
|Best et aL|p005| l we know that locally, despite the wide distribu¬ 
tions in host stellar mass of both HERGs and LERGs, LERGs have 
a strong preference to be hosted by galaxies with higher masses, 
while HERGs are hosted by galaxies with a lower median stellar 
mass but with a broader distribution. Assuming this still holds at 
higher redshifts, it implies that this strongly evolving population 
of lower mass Radio-Loud AGN are HERGs. We also expect the 
highest mass, most powerful sources to peak in space density at 
higher redshifts (2 ~ 2 — 3 |Rigby et al.|201 f| ) in line with the cos¬ 
mic downsizing picture where the most massive black holes have 
formed by 2 ~ 4. Indeed, results from many of the earlier radio sur¬ 


veys show that the most powerful (P^ 4 riH 7 ~ WHz ra¬ 


dio galaxies at 2 > 1 (e.g.jEales et al.|1997||jarvis et al.|2001 |Sey-| 


|mour et al.|2007[|Femandes et al. 12015 1 are predominantly HERGs 
hosted by the most massive galaxies. 

Indications from studies out to 2 < 1 show that the HERGs are 
indeed evolving more strongly with redshift than the LERG popu¬ 
lation ( [Best et al.|2014[ l such that within the redshift ranges of this 
study the radio-AGN population should be dominated by HERGs, 
opposite to that within the local universe. The mass dependence that 
we observe in this study supports this idea. Moreover, the evolution 
in the optical quasar luminosity function (i.e. that corresponding 
to radio-quiet cold mode accretion; |Hasinger et al.|2005||Hopkins| 
|et al.|2007{|Croom et al.|200^ is comparable to the kind of increase 
we observe for the low mass galaxies. 

Furthermore, the slope of the radio-loud fraction (cf. Fig. 
we observe at the higher redshifts, fuL oc is more con¬ 

sistent with the slope of the radio-loud fraction found in the local 
universe for HERGs only ^Janssen et al.|201^ . On the other hand, 
in the local sample, our derived slope of fuL oc matches 

that which [jMissen et ar] ( |2012^ found for only LERGs and matches 
the theoretical value for the accretion of hot gas from a halo (|Best| 
|et al.|2006) . In this case we know that the dominant population of 
all local radio sources is that of the LERGs. This lends support to 
the idea that there is an increase in the prevalence of HERG ac¬ 
tivity or cold mode accretion and that this mode is becoming the 
dominant population out to redshifts of 0.5 < 2 ^ 2. 

Lastly, as [Heckman & B^ ( |2014^ state, the crucial require¬ 
ment for cold mode accretion is the abundant central supply of cold 
dense gas. And as reflected in the increase in cosmic star formation 
rate density which has increased tenfold out to 2 ~ 2 and threefold 
out to 2 ~ 0.5 (e.g. jSobral et al.|2013[|Madau & Dickinson|2014| 
and references therein), there is significantly more cold gas fuelling 
star formation at these epochs. 

A closer inspection of Fig.[^reveals some interesting features, 
most particularly in highest mass bin, 11.25 < log(M*/M©) < 
12, in the highest redshift interval, 1 < 2 1.5. We suggest that 

while the radio-AGN sample as a whole at these redshifts is be¬ 
coming dominated by cold mode accretors, the highest mass, inter¬ 
mediate power ( 10 ^"* < f’]^ 4 GHz ^ 10 ^® WHz“'^) sources could 
still be hot mode/LERG sources because of their high mass. We 
observe a slight decrease in the number density of these specific 
sources consistent with the modelling results of jRigby et al.j where 
the ~ 10^^ WHz“^population peaks at 2 ~ 1 and is likely dom¬ 
inated by LERGs. This could be significant in the context of the 
models of LERG evolution presented byjBest et al.H2014|>. If these 


are indeed LERGs the decrease in number density does not match 
the order of magnitude decrease predicted by the preferred model 
of jBest et al.H2014} extrapolated out to 2 ~ 1.5 from fits to data out 
to 2 ~ 0.7. This model includes a time delay of ~ 1.5 Gyr between 
the formation of massive quiescent galaxies and when they are able 
to produce hot mode AGN, and at redshifts above ~ 1 the available 
population of host galaxies is declining so rapidly that such a delay 
in the onset of hot mode AGN activity would imply a drastic fall in 
the number densities of these sources above 2 ~ 1 . 


6 SUMMARY AND CONCLUSIONS 


In this paper we have used the SDSS value-added spectroscopic 
sample of radio-loud galaxies jBes^^^^c^mm 201 ^^and the 


VLA-COSMOS radio sample jSchinnerer et a . 2004| |2007[ l 


matched to a TTs-selected catalogue of the COSMOS/UltraVISTA 
field ( jMuzzin et al.|2013) to compile two samples of Radio-Loud 
AGN going out to 2 = 2. The samples are of sufficiently high radio 
power that they are dominated by RL AGN. Using these samples 
we have constructed radio luminosity functions for four host stellar 
mass bins between log(M*/M©) = 10.0 and log(M,/M 0 ) = 
12.0, in four redshift bins between 2 = 0.01 and 2 = 2. We 
have also investigated the radio-loud fraction as a function of stellar 
mass in these redshift bins. Together, we found the following: 


(i) Radio activity among galaxies of different mass increases 
differently towards higher redshifts. By considering the relative co¬ 
moving space density of radio-loud sources with respect to the lo¬ 
cal comoving space density, we showed that at 1.5 < 2 < 2 the 
space density of galaxies with stellar masses in the range 10.00 < 
log(M*/M 0 ) ^ 10.75 and 10.75 < log(M,/M 0 ) ^ 11.0 host¬ 
ing Radio-Loud AGN with ^ GHz > lO^'^WHz ^is respec¬ 
tively 45 ± 11 and 17.3 ± 2.6 times greater relative to the local 
space density while that of higher mass galaxies hosting Radio- 
Loud AGN remains the same. 

(ii) The fraction of galaxies which host Radio-Loud AGN with 

Pj 4 G]q 2 > 10 ^‘*WHz“^as a function of stellar mass shows a 
clear increase with redshift and a flattening with mass with the mass 
dependence evolving from /rl oc in the local sample to 

/rl oc Mb° at 1.5 < 2 < 2. 


We have argued that this increase in the prevalence of radio 
activity among galaxies of lower mass at higher redshifts is largely 
due to a rising contribution of AGN accreting in the radiative mode 
(HERGs). With this data we cannot yet conclusively show the evo¬ 
lution of the different accretion modes as a function of their host 
stellar mass because we lack the information on the excitation state 
of these sources at higher redshifts. However, future work combin¬ 
ing this and other radio-optical samples will allow more detailed 
studies of optical hosts of the high-redshift (1 < -z < 2 ) population 
of radio-AGN. 
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This study uses a Kj,-selected catalogue of the COS¬ 
MOS/Ultra VISTA field from jMuzzin et al.| ( |2013^ . The catalogue 
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contains PSF-matched photometry in 30 photometric bands cover¬ 
ing the wavelength range 0.15p.m —>■ 24p.m and includes the avail¬ 
able GALEX ([Martin et al.||2005|, CFHT/Subaru jCapak et al. 
2007l l, Ultra VIS' ^McCracken et al.|2012> , S-COSMOS [Sanders 


et al. 2007| l, and zCOSMOS ([Lilly et al.|2009|( datasets. This work 


uses the catalogue compiled by IBest & Heckman 

j2012|, which 

combines data from the VASC dBnnchmani^^ 

.l[2004[l of the 


SDSS DR7 ([Abazajianetal.[2009[l with NVSS ([Condon etal.|1998J 
and FIRST Isecker et al.[1995t . 
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